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ABSTRACT 

The transiting hot Jupiter planet Qatar-1 b was presented to exhibit variations in transit times that could be of perturbative nature. 
A hot Jupiter with a planetary companion on a nearby orbit would constitute an unprecedented planetary configuration, important 
for theories of formation and evolution of planetary systems. We performed a photometric follow-up campaign to confirm or refute 
transit timing variations. We extend the baseline of transit observations by acquiring 18 new transit light curves acquired with 
0.6-2.0 m telescopes. These photometric time series, together with data available in the literature, were analyzed in a homogenous 
way to derive reliable transit parameters and their uncertainties. We show that the dataset of transit times is consistent with a linear 
ephemeris leaving no hint for any periodic variations with a range of 1 min. We find no compelling evidence for the existence 
of a close-in planetary companion to Qatar-1 b. This finding is in line with a paradigm that hot Jupiters are not components of 
compact multi-planetary systems. Based on dynamical simulations, we place tighter constraints on a mass of any fictitious nearby 
planet in the system. Eurthermore, new transit light curves allowed us to redetermine system parameters with the precision better 
than that reported in previous studies. Our values generally agree with previous determinations. 


1. Introduction 


The transiting plane t Qatar-1 b (lAlsubai et al.ll201 ll) is the first planet ary body discovered by the Qatar Exoplanet Survey 


jAlsubai et al. 2013 ). It is a typical hot Jupiter on a 1.42 d circu lar orbit. 


Alsubai et al 


Coyinoetal 


(1201 11) give a mass of 1.09 ± 0.08 Mjup 


(120131) show that the planet has the greater mass 


(Jupiter’s mass), but new radial velocity (RV) measurements by 
of 1.33 ± 0.05 Mjup. The published studies provide consistent determinations of a planetary radius which is greater by 16-18% 
than Rjup (Jupiter’s radius). The planet produces transits which are nearly 1.5 h long and 2.3% deep in flux. The host star, GSC 
4240-470 {V - 12.84 mag), is a metal-rich d warf of spectral typ e K3. It exhibits a moderate chromospheric activity, which 
could be partly induced by the close-in planet ( Covino et al.ll2013 ). The RV observations throu ghout a transit show a signature 
of the Rositter-McLaughlin effect that is consistent with a sky-projected obliquity close to zero ( Covino et al.ll2013 ). 


More recently, 


von Essen et al. 


(120131) have found indications for long-term variations in transit timing residuals from a linear 


ephemeris. The periodicity of this signal is postulated to be close to 190 or 390 d with a peak-to-peak amplitude of about 2 min. 


*Partly based on (1) data collected with telescopes at the Rozhen National Astronomical Observatory and (2) observations obtained with 
telescopes of the University Observatory Jena, which is operated by the Astrophysical Institute of the Friedrich-Schiller-University. 






























Table 1. Details on instruments taking part in the campaign. 


# 

Telescope 

Observatory 

Diameter (m) 

CCD detector 

Size of matrix and field of view 

Mr 

1 

Ritchey-Chretien-Coude Telescope 

National Astronomical Observatory Rozhen, Bulgaria 

2.0 

Roper Scientific VersArray 1300B 
1340 X 1300 pixels, 518 x 516 

1 

2 

Ritchey-Chretien Telescope 

Sierra Nevada Observatory (OSN), Spain 

1.5 

Roper Scientific VersArray 2048B 
2048 X 2048 pixels, 7192 x 7192 

4 

3 

Trebur 1-meter Telescope 

Michael Adrian Observatory, Trebur, Germany 

1.2 

SBIG STL-6303 

3072 X 2048 pixels, 1010 x 617 

3 

4 

Schmidt Teleskop Kamera (Mugrauer & Berthold 2010) 
University Observatory Jena, Germany 

0.9 

E2V CCD-42-40 

2048 X 2048 pixels, 5218 x 5218 

3 

5 

Cassegrain Telescope 

Nicolaus Copernicus University, Poland 

0.6 

SBIG STL-1001 

1024 X 1024 pixels, 1118 x 1118 

6 

6 

Cassegrain Telescope 

Volkssternwarte Kirchheim, Germany 

0.6 

SBIG STL-6303 

3072 X 2048 pixels, 5214 x 3419 

1 


Notes. Diameter is the size of a telescope’s main mirror. N„ is the total number of complete or partial transit light curves, acquired with the 
given instrument. 


Table 2. Details on individual light curves reported in this paper. 


Date UT (epoch) 

Telescope 

Filter 

r 

pnr 

A 


Sky conditions 

2011 Aug 01 (181) 

Volk.Kirch. 0.6 m 

none 

0.87 

3.68 

1.05 ^ 1.03 ^ 1.05 

Clear 

2012 Dec 13 (533) 

Torun 0.6 m 

none 

1.06 

2.43 

1.12 ^ 

1.55 

Clear 

2013 Sep 03 (719) 

Torun 0.6 m 

none 

1.51 

2.51 

1.02 ^ 

1.14 

Clear 

2013 Sep 30 (738) 

Torun 0.6 m 

none 

1.24 

2.09 

1.03 ^ 

1.22 

Photometric 

2013 Nov 26 (778) 

Trebur 1.2 m 

none 

1.39 

1.18 

1.07 ^ 

1.14 

Photometric 

2013 Dec 20 (795) 

Jena 0.9 m 

Rb 

1.12 

2.67 

1.23 ^ 

1.81 

Clear 

2014 Mar 25 (862) 

Trebur 1.2 m 

none 

0.40 

4.36 

2.10 ^ 

1.43 

Passing clouds 

2014 Apr 02 (867) 

Rozhen 2.0 m 

Rc 

1.22 

1.47 

1.60 ^ 

1.22 

Photometric 


Torun 0.6 m 

none 

1.82 

2.36 

1.41 ^ 

1.14 

Clear 

2014 May 16(898) 

OSN 1.5 m 

Rc 

1.09 

2.21 

1.48 ^ 

1.14 

Mostly clear, some high clouds 

2014 May 26 (905) 

OSN 1.5 m 

Rc 

1.33 

1.39 

1.62 ^ 

1.16 

Mostly photometric 

2014 Jun 22 (924) 

OSN 1.5 m 

Rc 

1.09 

1.61 

1.37 ^ 

1.13 

Mostly photometric, occasionally 
passing high clouds before ingress 

2014Jul01 (931) 

Trebur 1.2 m 

Rb 

0.88 

2.64 

1.17 ^ 

1.03 

Photometric, first half of series 
partly obscured by the dome slit 

2014 Sep 10(981) 

OSN 1.5 m 

Rc 

1.33 

1.29 

1.14^ 1.13 ^ 1.37 

Photometric 

2014 Sep 18(986) 

Jena 0.9 m 

Rb 

1.16 

2.25 

1.19 ^ 

1.66 

Photometric 

2014 Sep 27 (993) 

Jena 0.9 m 

Rb 

1.14 

2.02 

1.12 ^ 

1.54 

Photometric 


Torun 0.6 m 

Rc 

0.95 

3.74 

1.14^ 

1.53 

Clear, high humidity 

2014 Oct 27(1014) 

Torun 0.6 m 

Rc 

0.86 

3.45 

1.05 ^ 

1.21 

High humidity, occasionally passing 
high clouds 


Notes. Date UT is given for mid-transit time. Epoch is the transit number from the initial ephemeris given in lAlsubai et alj l l2011h . Rc and 
Rb denote Cousins and Bessell f^- band filters, respecti vely. E is the median number of exposures per minute, pnr is the photometric scatter in 
millimag per minute as defined bv iFulton et al.l l l201 ih . X shows the airmass change during transit observations. 


These variations were interpreted as a result of the exchange of energy and angular mo mentum via gravitational interaction 
with an unseen, close-in planetary companion ( Holman & Murra^boOS lAgol et al.ll20Q5ll . Such planetary configuration would 
be unprecedented because hot Jupiters appear to be alone or accompanied by other planetary bodies on wide orbits (e.g., 


Steffen et al 


2012h . On the other hand, studies of candidates for hot planets from a Kepler s ample show that s ome fraction of 


these objects reveal periodic variations in transit timing that could be dynamically induced (ISzabo et al.ll2013ll . There may be 
only a small fraction of hot Jupiters which are a part of compact multiple-planet systems. If such systems are confirmed, their 
properties will place constraints on planet formation and migration theories. In some circumstances, migration models, based 
on tidal interactions between a giant planet and a circumstellar gas disk, predict the presence of ne ighboring low-mass planets, 
close to or trapped in interior or exterior mean-motion resonances (MMRs) with hot Jupiters (e.g., 


Podlewska & Szuszkiewicz 


20091). Such resonant or near-resonant companions would induce orbital perturbations producing detectable variations in transit 


times. The analysis of multiple-planet candidate systems shows that most planetary systems ar e not resonant, but the distribution 
of planet period ratios shows clumping just wide of low-order resonances, such as 2:1 or 3:2 ( Fabrvckv et al I2OI4I) . The origin 
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Fig. 1. New transit light curves acquired for Qatar- 1 b, sorted by observation dates. The continuous lines show the best-fitting model, the 
residuals are plotted below each light curve. 

of this pattern is unclear. The question whether hot Jupiters are accompanied by planets on nearby orbits or not, and what is the 
fraction of such planetary configurations still remains open. 

In this context, Qatar-1 b with its preliminary detection of the transit time variation (TTV) signal could become a prototype 
of a new group of multiple-planet packed systems with hot Jupiters. As the system deserves further investigation, we organized 
an observing campaign to confirm or refute perturbations in the orbital motion of Qatar-1 b by the method of transit timing. 


2. Observations 

We acquired 18 new light curves for 16 transit events with six 0.6-2.0 m telescopes located across Europe (in Bulgaria, Poland, 
Germany, and Spain). The characteristics of individual instruments and detectors are given in Table[T] Two transits were observed 
simultaneously from two distinct sites. Most of the observations were performed in R filter, in which the instrumental setups 
were found to be the most effective. We also acquired light curves in white light, i.e. without any filter, to increase the signal- 
to-noise ratio and obtain more precise data for timing purposes. Telescopes were moderately defocused to observe the Qatar-1 
star with exposure times between 10 and 60 s, depending on mirror diameter, filter, and observing conditions. For a given light 
curve, the exposure time was kept fixed throughout a run. The timestamps were saved in coordinated universal time (UTC), 
provided by a GPS system or the network time protocol software. Details on individual runs are presented in Table |2] 










































Table 3. Redetermined system parameters, together with values from the literature. 


Parameter 

Orbital inclination, (deg) 
Scaled semi-major axis, 

Planet to star radii ratio, RbjR* 
Transit parameter, b = ^ cos ib 
Transit total duration, T 14 (min) 
Semi-major axis, a\, (AU) 
Planetary mass*, Mb (Mjup) 
Planetary radius, 7?b {Rjup) 
Planetary density, pb (pjup) 
Planetary gravity*, gb (m s“^) 
Stellar mass, M, (Mq) 

Stellar radius, R, (Rq) 

Stellar density, p, (pg) 

Stellar gravity‘', logg, (cgs) 
Stellar luminosity, L, (Lg) 

Age, (Gyr) 


This paper 
84.26^;' 

0 14 ^ 01+000076 
'^■^^-'^^-0.00078 

0.63 + 0.02 
98.5 + 1.7 

A A99QQ+0.00069 

0-87io! 

25.65 + 0.71 

0.803^00^2 

0.782!p| 
1.680!|:0| 
4.556 
0.318 
9.8 


0.016 

0.051 

■'- 5.8 


Alsubai etal. (2011)“ 


83.47-!^ 


0.1453 + 0.0016 

0 696^0021 
u.o;?o_qq24 


0 02343 + 0.00026 
u.uz,j‘tj_QQon25 

1 090^008^ 
i.Ub'U gQgj 

1.164 + 0.045 

0 690^0098 
u.<Jb'U_QQ84 
18 4+20 

0.85 + 0.03 
0.823 + 0.025 
1.52 + 0.12 
4.536 + 0.024 


>4 


Co vino et al. (2013) von Essen et al. (2013) 


83.82 + 0.25 
6.25 + 0.10 
0.1513 + 0.0008 
0.675 + 0.016 
97.6 + 1.4 
0.0234 + 0.0012 
1.33 + 0.05 
1.18 + 0.09 
0.80 + 0.20 
23.6:|'4 
0.85 + 0'.03 
0.80 + 0.05 
1.62 + 0.08 
4.55 + 0.10 


84.52 + 0.24 
6.42 + 0.10 
0.1435 + 0.0008 


:4.5 


Notes. Values for a solution with a circular orbit were taken. ® In calculations, the value of the RV semi-amplitude was taken from Covino 
et al. (2013). Based on parameters derived from transit photometry. 


Photometric time series were obtained with the differential aperture photometry applied to science images after bias and dark 
subtraction followed by flat-field correction with sky flats. The radius of the photometric aperture was optimized for each dataset 
to achieve the smallest scatter in the out-of-transit light curves. Magnitudes were determined against available comparison stars. 
Each light curve was de-trended by fitting a parabola along with a trial transit model, then converted to fluxes and normalized 
to unity o utside of the transit event. The timestamps were converted to barycentric Julian dates in barycentric dynamical time 

' ilbninl) . 


(BJDtdb, 


Eastman et al.ll 


3. Results 


3.1. System parameters 

The Transit Analysis Package (TAP, 


Gazak et al 


1l2ni2ll 


was used to model the s et of 1 8 new transit light curves, supplemented 


Covino et al. 


(120131) provide 5 high-quality photometric time 


with 31 complete or partial light curves available in the literature. 
series acquired betwee n May 2011 and September 2012 with the 1.82 m Asiago (Italy) and 1.23 m Calar Alto (Spain) telescopes. 


von Essen et al 


(120 13 I) make available 26 light curves. They were obtained from March 2011 to August 2012 with the 1.2 m 


Oskar-Liihning telescope at Hamburg Observatory (Germany) and the 0.6 m Planet Transit Study Telescope at the Mallorca Ob¬ 
servatory (Spain). The TAP code uses the Markov Chain Monte Carlo (MCMC) method with the Metropol i s-Has tings algorithm 
and a Gibbs sampler to find the best-fit transit model, based on the analytical approach of lMandel & Ago] (2002). As the photo¬ 
metric time series may be afflicted by time-correlated noise, the wavelet-based alg orithm of Carter & WinnI (l2009h is employed 
to robustly estimate parameter uncertainties. A quadratic limb darkening (ED) law jKonalll950l) is implemented to approximate 
the flux distribution across the stel lar disk. The values of li near and quadratic ED coefficients, mi and M 2 , res pectively, were 


linear ly interpolated from tables of IClaret & BloemenI (l201 ll) using a dedicated tool of the EXOEAST package dEastman et al 


20131). The stellar parameters, such as effective tempe rature r^fr - 4910+ 100 K, surface gravity logg = 4.55+0.1 (in cgs units), 
and metallicity [Ee/H] = 0.2 + 0.1, were taken from ICo vino et al.l (120131) . The efficient maximum of the energy distribution in 


white light was found to coincide with the R band, so J?-band ED coefficients were adopted for white light curves. 

In a final iteration of the fitting procedure, the orbital inclination fb, the semimajor-axis scaled by stellar radius a\,fRt, and 
the planetary to stellar radii ratio R\,/Rt were linked for all light curves and fitted simultaneously. The mid-transit times were 
determined independently for each epoch. The orbital period Pb was kept fixed at a value refined in Sect. 13.21 The ED coefficients 
were allowed to vary around the theoretical values under the Gaussian penalty of cr = 0.05 to acco unt for uncertainties in stellar 
parameters and possible systematic errors of the theoretical predictions (e.g., iMiiller et al.ll2013h . The orbit of the planet was 
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Fig. 2. Modified Hertzsprung-Russel diagram with Qatar- 1 marked as a central dot. The PARSEC isochrones of [Fe/H] = 0.20 and ages 
between 3 and 13 Gyr with a step of 2 Gyr are sketched with dashed lines. The best-fitting isochrone is drawn with a continuous line. 


assumed to be circular (ICovino et al.ll2013ll . Although the individual light curves are de-trended, the TAP code was allowed to 
take possible linear trends into account and include their uncertainties in the total error budget of the fit. We used ten MCMC 
chains, each containing 10® steps. The first 10% of the results were discarded from each chain to minimize the influence of 
the initial values of the parameters. The best-fitting parameters, determined as the median values of marginalized posteriori 
probability distributions, and 1 cr uncertainty estimates, defined by a range of 68.3% values of the distributions, are given in 
Table[3] The new light curves are plotted in Fig. [T] together with the transit model and residuals. 

The parameters derived from the transit model, combined with literature quantities determined from spectroscopic data, 
allowed us to redetermine the stellar mass M,, luminosity L,, and system age. The mean stellar density p, was calculated from 
transit parameters, independent of theoretical stellar models, using the formula 


Stt / flb 


GPl\R, 


( 1 ) 


where G is the universal gravitational constant. Then, the host star was placed at a modified He rtzsprung-Russel dia gram 


(Fig. |2ll displaying the p, versus Teff plane, together with PARSEC isochrones in version 1.2S (IBressan et al.ll2Q12h . The 


stellar parameters and their uncertainties were derived by interpolating the sets of isochrones with [Fe/H] ranging from 0.1 to 
0.3. 

The parameters and Rb/R, were checked for long timescale trends or periodic variability (abIRt was not considered 
because it is correlated with ib). The fitting procedure was repeated for a subset of 37 complete light curves, with ib and Rb/R* 
allowed to be determined independently for each dataset. The parameter Ub/Rt was allowed to vary around a best-fit value 
under the Gaussian penalty, as given in Table |3] Individual determinations are plotted in Fig. |3] A model assuming a constant 
ib has a reduced chi-squared of 1.15 and a p-value of 0.25. The parameter Rb/R, is also constant with 
corresponds to a p-value of 0.15. For both parameters, there is no reason to reject a null hypoth esis assuming constancy. In 


addition, both parameters were searched for periodic variations with the Fomb-Scargle algorithm (IFomb 


1976 


Scarglelll982h 


and no significant signal was found. The lack of any variations justifies linking these parameters in the final model. 


3.2. Transit timing 

New mid-transit times, together with redetermined values from the literature data, allowed us to refine the transit ephemeris. As 
a result of a least-square linear fit, in which timing uncertainties were taken as weights, we derived 

Pb = 1.42002406 + 0.00000021 d (2) 

and the mid-transit time at cycle zero 


To = 2455518.41117 +0.00014 BJDtdb ■ 


(3) 
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Fig. 3. Distribution of and in a function of time for a trial fit, in which both parameters were allowed to vary for individual complete 

transit light curves. Open and filed symbols represent literature and our new datasets, respectively. Continuous lines mark weighted mean 
values. Grayed areas between dashed lines show uncertainties at the confidence level of 95.5%, i.e. 2cr, where cr is the weighted standard 
deviation. 
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Fig. 4. Residuals of transit times from the refined linear ephemeris. Open symbols mark literature data. Filled dots denote new transits reported 
in this paper. T he mid-transit time at e poch 181 (2011 Aug 01 = BJD 2455775.4), marked with a half-filled symbol, is determined using a 
light curve from ivon Essen et"^ ll2013h and a light curve reported in this paper. The grayed area between dashed lines shows the propagation 
of uncertainties of the ephemeris at a 95.5% confidence level. 


The value was found to be equal to 1.07, clearly showing that there is no hint for any deviation from the constant orbital 

period. The procedure was repeated for a limited sample of 17 mid-transit times with uncertainties smaller than 40 s, hence the 
most reliable. Again, the linear ephemeris was found to reproduce observed transits perfectly with;^^^^^ = 0.91. This shows that 
any possible TTV signal is not masked by a scatter of data points of lower quality and hence lower timing precision. Individual 
transit times and residuals are collected in Table|4l and residuals from the ephemeris are plotted in Fig.|4] 

The Lomb-Scargle periodogram (Fig. |5j reveals no significant power at periods of 187 + 17 and 386 + 54 d, reported by 
( 12013h . A bootstrap resampling method, based on 10^ trials of the randomly permuted timing residuals at the 


von Essen et al 


original observing epochs, was used to determine a false alarm probability (FAP) of 28% for the strongest peak. 


3.3. Constraints on an additional planet 

The TTV method is particularly sensitive to perturbers close to MMRs. The non-detection of periodic variations in transit timing 
allows such planetary configurations to be ruled out. In remaining configurations, planetary companions may be detected easier 
in the domain of precise Doppler measurements. Combining these two methods together places constraints on an upper mass of 
any hypothetical second planet in the system as a function of its semi-major axis. 
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Table 4. New and redetermined mid-transit times. 


Epoch 

91 

131 

136 

162 

169 

181 

198 

200 

224 

329 

367 

369 

374 

381 

393 

407 

417 

419 

426 

431 

438 

450 

455 

467 

469 

481 

502 

533 

719 

738 

778 

795 

862 

867 

898 

905 

924 

931 

981 

986 

993 

1014 


Emid (d) 


56405338”^“ 

5647 . 63267 +°°“ 5 ’^ 


5704.43426 
5711.53450 
5748.45576!' 
5758.39484! 
5775.43517! 
5799.57630! 
5802.41602! 
5836.49672! 
5985.6005!° oojj 


-o.c 


6039.56043+° 
6042.40053 
6049.50010! 
6059.4393!°- 
6076.4792 
6096.36008, 
6110.56029! 
6113.4004!°- 
6123.34169! 
6130.44153! 
6140.38094! 
6157.42152+ 
6164.52187 
6181.56264! 
6184.4027!°- 
6201.44229! 
6231.26318+ 
6275.28493 
6539.40816 
6566.38885+ 
6623.18947+ 
6647.33094! 
6742.4720!°- 
6749.57192! 
6793.59237! 
6803.53269! 
6830.51327+ 
6840.45424 
6911.45488+ 
6918.55468+ 
6928.49510! 
6958 3143 ^^- 


00075 


076 


)0035 

)0052 

II? 


>0043 

>0064 


O-C (d) 

-hO.0006 

-0.00069 

-0.00005 

-hO.00006 

-hO.00070 

-0.00039 

-0.00036 

-hO.00037 

-hO.00004 

-hO.00017 

-hO.0014 

-hO.00043 

-hO.00049 

-0.00007 

- 0.0011 

-0.0015 

-0.00087 

-0.00091 

-0.0009 

-hO.00028 

- 0.00001 

-0.00077 

-0.00047 

-0.00024 

-hO.00024 

-hO.0003 

-0.00045 

-0.00006 

-hO.00094 

-0.00030 

-0.00007 

-0.00042 

-hO.00065 

-hO.00014 

- 0.00011 

-0.00040 

-0.00025 

- 0.00012 

+0.00068 

+ 0.00012 

- 0.00020 

+0.00004 

-0.0013 


Data source 

1 

1 

1 

1,2 

1 

1 

1,3 

2 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 

1,2 

1 

1 

2 

2 

1 

1 

1 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 


Notes. Epoch is the transit number from the initial ephemeris given in lAlsubai et alj (1201 Ih . is the mid-transit time in B JDtqr minus 
24500 00 d. O-C is the value of timing residuals from the refined transit ephemeris. Data source: tll lvon Essen et al]ll2013h : (2l lCovino et alJ 
( 120131) : (3) this paper. 


We followed the method which we applied to WASP-3 and WASP-1 systems ( Macieiewski et al.l2013 , 2014 ). The Bulirsch- 
Stoer integrator, implemented in the Mercury 6 code ( Chambersll999 ). was used to predict deviations from a Keplerian motion 
of Qatar-1 b caused by a fictitious perturbing planet. The system was assumed to be coplanar, with both orbits initially circular. 
The mass of the perturber was set at 0.5, 1, 5, 10, 50, 100, and 500 Msarth (Earth masses), and the initial semi-major axis varied 
between 0.023 and 0.107 AU (corresponding to the orbital period ratio between 1 and 10) with a step of 2 x 10^° AU. Only 
configurations with outer perturber were considered because of the short orbital period of Qatar-1 b. The initial orbital longitude 
of the transiting planet was set at a value calculated for cycle zero, and the initial longitude of the fictitious perturber was shifted 
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Fig. 5. Lomb-Scargle periodogram for timing residuals plotted in Fig. l^ Dashed and dotted h orizontal lines show the 5% and 1% FAP levels, 
respectively. Arrows indicate positions of the TTV signals claimed bv ivon Essen et 



Fig. 6. Upper mass limit for a fictitious outer planet in the planetary system around Qatar-1, based on timing and RV datasets, as a function 
of the orbital period of that planet, Pp. Configurations that are below the detection threshold of both techniques are placed in the grayed space 
of parameters, bounded by a continuous line. The dotted and dashed lines show detection limits for timing and RV techniques, respectively. 
Configurations with Pp/Pb < 1-5 were found to be unstable in a short timescale of the simulations. 


by 180°. The calculations covered a time span of 1500 days in which transits of Qatar-1 b were observed. The residuals from 
a linear ephemeris for synthetic datasets were compared with the observed nns of transit timing, and a TTV upper mass of a 
fictitious planet was calculated. 

Re-analysis of the RV dataset from 


Covino et al 


(120131) . performed with the Systemic 2.16 (IMeschiari et al.ll2009l) . gives 
a circular single-planet orbit with nns^y = 11.1 m s '. The RV measurements were acquired with the High Accuracy Radial 
velocity Planet Searcher-North (HARPS-N) spectrograph installed at the 3.58 m Tele scopio Nazionale Ga lileo (TNG) on La 


Palma (Spain). The RV uncertainties are between 4.4 and 20 m s *. The dataset from 


Alsubai et al 


(1201 ih . acquired with the 


Tillinghast Reflector Echelle Spectrograph (TRES) coupled with the 1.5-m telescope at the Ered L. Whipple Observatory (USA), 
was omitted because of the much lower precision of the RV measurements (22 -61ms'*). The value of nnsj^v was transformed 
into the RV mass limit of the hypothetical second planet on a circular orbit as a function of its semi-major axis. 

Einally, the TTV and RV criteria were combined together to find planetary configurations with a hypothetical second planet 
below the detection threshold. Results are illustrated in Eig. |6] The TTV technique is more sensitive for configurations within 
the 2:1 resonance. We can eliminate perturbers with masses between 10 and 20 Meanh out of MMRs, and probe for a sub-Earth 
mass regime in 2:1 and 5:3 MMRs. For orbital periods longer than about 2.2 P^, the Doppler dataset provides tighter constraints 
except a 3:1 MMR, in which a perturber with a mass down to about 3 Mnarth can be eliminated thanks to the TTV method. 


4. Concluding discussion 


Despite the conservative approach, our light curve modeling provides system parameters with the smallest uncertainties pub¬ 
lished so far. This is the result of the homogenous analysis of the rich set of 18 new and 3 1 literature light curve s. The or¬ 
bital inclination /b = 84?26''''’°'^ 
fers by 1.1 and 1.4 cr from results of 


agre es perfectly w i th the valu e of 84?52 + 0?24 re ported by 


von Essen et al 


Covino et al 


(120131) and 


Alsubai et al 


flb/R* = 6.319 


+ 0.070 

- 0.068 


is between 6.25 + 0.10 and 6.42 + 0.10 brought by 


_ _ (120131) . and dif- 

( 201 1[), respe ctiv ely. The scaled semi-m ajor axis 


Covino et al. 


(120131) and 


von Essen et al. 


tively. The planet to star radii ratio Rb/R* - 0.1459NQQyQ7g is consistent with 0.1453 + 0.0016 derived by 


(l2013l). 


Alsubai et al 


respec - 

( 12011 1 ) . 
















































































and differs by 1.5 and 3.5 cr from values reported by 
physical properties derived in this study confirm results of 


von Essen et al 


( 2013 )and 


Covino et al. 


Covino et al 


Alsubai et al 


(120 13h . respectively. Planetary 
(120 13h . and deviate slightly from values reported by 
(1201 Ih due to the planetary mass underestimated in that study. Stellar parameters were found to agree with values 
determined in previous studies. The system’s age is 10 Gyr with the lower uncertainty of a; 6 Gyr and the upper one limited 
by the age of the Universe. The host star is in the middle of its lifetime on the main sequence that is estimated to be 18 Gy for 
a 0.8 Mq star. 

Spectral obser vations of the Ca II H and K lines show that the host star exhibits a moderate chromospheric activity with 
logfJJ^j, = -4.60 (ICovino et al.ll2013l) . However, we find no statistically significant variation in planet to star radii ratio that 
could be c aused by the evolvi ng distribution of dark spots or bright faculae outside the projected path of the planet on the stellar 


Czesla et al 


20091). We also identified no features in the transit light curves which could be attributed to starspot 


disk (e.g., 
occultations by the planetary disk. 

Transit timing observations clearly show that the orbital motion of Qatar-1 b is not perturbed by any body which could 
produce periodic deviations with a ra nge greater than 1 minu te. We find no evidence which confirms the TTV signal with the 


periodicity and amplitude reported by 


von Essen et al 


(120131) . Combining the transit timing and RV datasets, we can rule out 


three proposed scenarios with perturbers in 2:1, 5:2, and 3:1 resonances even for circular orbits (in general, eccentric orbits 
make perturbations more pronounced). The lack of the TTV signal also makes a proposed massive perturber on a wide orbit 
unlikely. We note, however, that further RV measurements, acquired in a course of a few hundred days, would unequivocally 
reject this scenario. 

We conclude that Qatar-1 b has no detectable planetary companions on nearby orbits. This finding is in line with results based 
on a sample of hot Jupiters obse rved with the space -borne facilities or ground based telescopes. As the formation of hot Jupiters 
is not yet fully understood (e.g., Steffen et al.ll2012 ). the loneliness of Qatar-1 b speaks in favor of i nward migration theories of 


massive outer planets thro ugh planet-planet scattering caused by mutual dynamical perturbations (IWeidenschilling & Marzari 
1 9961: iRasio&Fordll 19961) . 
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